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We study how photon absorption losses degrade the bipartite entanglement of entangled states of 
light. We consider two questions: (i) what state contains the smallest average number of photons 
given a fixed amount of entanglement? and (ii) what state is the most robust against photon 
absorption? We explain why the two-mode squeezed state is the answer to the first question but 
not quite to the second question. 



I. INTRODUCTION 

For long-distance quantum communication the most 
important decoherence mechanism is due to photon ab- 
sorption. It would be nice to know what type of entangled 
states of light is most robust, i.e. preserves its entangle- 
ment best, against photon absorption losses. Intuitively 
one might think that states containing fewer photons will 
be more robust. After all, both loss rates and decoher- 
ence rates tend to increase with the number of photons. 
For instance, a cavity filled with exactly N photons loses 
photons at a rate proportional to N . Similarly, it is well 
known that, in general, the more macroscopic a state is 
the more sensitive it is to noise. For example, the deco- 
herence rate of Schrodinger cat states of the (unnormal- 
ized) form 

\a)±\-a), (1) 

with I a) a coherent state, increases with the amplitude 
\a\. This phenomenon was demonstrated experimentally 
in the context of microwave cavity-QED Q and for mo- 
tional states of ions in an ion trap (2|]. Also entangled 
coherent states of the (unnormalized) form 

\a)\a)-\-a)\-a), (2) 

decohere faster with increasing \a\. This fact was dis- 
cussed in several theoretical papers: Ref. uses an op- 
erational measure and shows that the teleportation fi- 
delity decreases monotonically with |a|, and another pa- 
per [Jl demonstrates that the entanglement of formation 
decreases monotomically with \a\ as well. 

If the intuition that a smaller number of photons makes 
a state more robust were true, then the most robust en- 
tangled state of light would be the two-mode squeezed 
state. That state contains the smallest average number 
of photons given a fixed amount of entanglement, as we 
will show in Section II. We will see, however, that the 
intuitive answer is not quite correct and we will explain 
why this is so in Section III. The calculations in Section 
III proceed as follows: we consider large classes of pure 



states (both Gaussian and non-Gaussian states) with a 
fixed amount of entanglement. Then we calculate the 
entanglement that remains after those states have deco- 
hered by being subjected to photon absorption. That 
remaining entanglement is considered to be a measure of 
the robustness of the entangled state against photon ab- 
sorption. The robustness is investigated as a function of 
the average number of photons in the initial pure state. 
We use two definitions of entanglement, the entanglement 
of formation and the negativity. In the former case we 
restrict our attention to states for which the entangle- 
ment of formation can be calculated analytically, in the 
latter case there is no restriction on the types of states 
considered. 

We end this Introduction by noting that the entangle- 
ment of a state of the electromagnetic field is completely 
determined by its expansion in Fock (photon number) 
states 5] . Consequently, the effects of photon absorption 
losses on the entanglement of a state are likewise com- 
pletely determined by the form of the state in the Fock 
basis. On the other hand, the effects of other decoher- 
ence mechanisms arising from, say, polarization diffusion 
or phase diffusion, depend on the physical implementa- 
tion of the state. For example, if one uses two modes that 
are different only in their polarization degree of freedom, 
then polarization diffusion mixes the two modes. But 
if the two modes are spatially separated then no such 
mixing occurs. The type of decoherence discussed in the 
present paper is thus the most universal type of deco- 
herence of entangled states of light, in the sense that it 
is independent of the precise physical implementation of 
the state. 



II. PHOTONS AND ENTANGLEMENT 

Let us consider a bipartite system consisting of two 
spatially separated electromagnetic field modes A and B. 
The entanglement in a pure state on A and B is given by 
the von Neumann entropy of the partial density matrix 
of either system A ot: B, where the pure state can be 



2 



written in the photon-number basis as 

l^')4.B= > anm\n) A\m) B ■ (3) 



}A,B = ^ 
n.m—Q 



The problem of finding the state with the smallest aver- 
age number of photons for a fixed amount of entangle- 
ment is mathematically equivalent to a well-known prob- 
lem from thermodynamics: finding the state with the 
largest entropy given a fixed energy, which is, almost ob- 
viously, the same as finding the state with the smallest 
energy given a fixed entropy. Indeed, energy is propor- 
tional to the average number of photons, and the entropy 
is just the entanglement. The answer to the thermody- 
namics question is of course the thermal state. That is, 
the state 



= (1-A)^A" 



(4) 



maximizes the entropy for fixed average energy, where 
A = exp(— ?iw/fcT) if nTwj is the energy of the state \n) 
and T the temperature. But that mixed state arises in- 
deed from a two-mode squeezed state of the form 



(5) 



for any 7 for which I7P = A. In particular, the entangle- 
ment in this state is 



E ={- + !) log2(n + 2) - nlog2 n - 1, 
in terms of the average number of photons 
- 2A 2|7p 



1-A 1 



I7|- 



(6) 



(7) 



For instance, for n = 1 one has more than one unit of en- 
tanglement, as = 31og2(3)/2 — 1 k, 1.3774. Conversely, 
it takes less than one photon (0.5876 of a photon to be 
more precise) to have one ebit of entanglement between 
two modes. That might seem surprising as one might 
have expected the state |0)|1) -f |1)|0), a state with one 
photon and one ebit of entanglement, to maximize the 
entanglement of one photon in two modes. 

If one allows an arbitrarily large number of modes for 
both parties, it is straightforward to get an arbitrar- 
ily large amount of bipartite entanglement. For exam- 
ple, consider states with two photons: with M modes 
for each party, the bipartite state 1 1000 ... 1000 .. .) -I- 
10100 . . . 0100 . . .) -I- 10010 . . . 0010 . . .) -I- . . . has exactly 
log2(M) ebits of entanglement. But the state that max- 
imizes the amount of entanglement in the case of M 
modes for both parties for two photons contains more 
than twice as much entanglement and is simply a tensor 
product of M two-mode squeezed states each with aver- 
age photon number 2/M, with a total entanglement of 
21og2(Af) + 1 + 0{l/M) for large M. 



III. ENTANGLEMENT OF DECOHERED 
STATES 

In this Section we will consider the most robust bipar- 
tite state of just two modes. This should reveal most 
properties of multi-mode entangled states under the as- 
sumption that all modes decohere independently and 
identically. In particular, we assume the following model 
for photon absorption losses: For any coherent state \a)A 
of a mode photon absorption is governed by the pro- 
cess 



|0)j 



\^a)AWl - m)i 



(8) 



where E denotes the environment, which is assumed to be 
unobservable and hence must be traced out. The param- 
eter rj gives the fraction of photons in mode A that sur- 
vives the photon absorption process, with ?? = 1 thus cor- 
responding to a noiseless channel. Since coherent states 
form an overcomplete set of states, the model © fully 
specifies decoherence due to photon absorption. We as- 
sume that all modes decohere in the same way, according 
to ©, with the same parameter for each mode and a 
different environment for each mode. 

A pure state of the form Q is turned into a mixture 
by the photon absorption process (jHJ. Calculating the 
entanglement of a general mixed state is no trivial prob- 
lem, as is well known. We follow two paths here in order 
to find the entanglement that remains after photon ab- 
sorption. 

In the first subsection we consider states for which we 
can analytically calculate the entanglement of formation 
for the mixed state that arises from photon absorption. 
This first of all includes all states in which each of the 
two modes A and B contain at most one photon. In that 
case the Hilbert space of each mode is two-dimensional, 
also after photon absorption, and we can apply the Woot- 
ters formula . Secondly, we can consider the two- mode 
squeezed state, as it is a Gaussian state and remains 
Gaussian after photon absorption, so that we can apply 
the results of ■ In the second subsection we consider ar- 
bitrary states and calculate (numerically) the negativity 
Here too, analytical calculations become tedious 
or impossible depending on the dimension of the matrix 
whose eigenvalues have to be calculated. 

In both subsections we consider states with a fixed 
amount of entanglement and wish to find the states that 
preserve their entanglement best under photon absorp- 
tion. 



A. Entanglement of formation 

We first consider states of the general form 



(9) 



where the parameter < p < 1/2 fully determines the 
entanglement of formation of the state 



E = -plog2ip) - (1 -_p)log2(l -p). 



(10) 
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FIG. 1: Entanglement of formation as a function of the av- 
erage number of photons (in both modes combined) for all 
bipartite states with at most one photon in each mode. 

and where 

\4')a = cosq;|0) + sina|l) 

\^p)b = cos/3|0) +sin/3|l) 

\4>^)a = — sin ajO) + cos 

\ip^)B = -sin/3|0) +cos/3|l) (11) 

form arbitrary orthonormal bases in the subspaces 
spanned by the zero- and one-photon states in modes 
A and B. We then fix p and 77 and plot the entanglement 
of all states of the form © after they have lost a frac- 
tion 1 — ?7 of their photons by varying the angles a and 
(3 independently. We plot the result as a function of the 
average number of photons n in the original undecohcred 
pure state 

n = p(cos^ a + cos^ /?) -f (1 - p)(sin^ a + sit? (3). (12) 

Results for specific parameters are plotted in Fig. 1. In 
particular, we choose 77 = 0.5 and fix p = 1/3, so that 
the entanglement of formation of the pure state is i? ~ 
0.9183. The optimum state is found to be 

V^|01)-v/T^|10), (13) 

which contains one photon and for which the entangle- 
ment left after decoherence is £' « 0.3236. Now clearly, 
this optimum is not obtained for the smallest average 
number of photons, which is 2/3 for the states con- 
sidered. That minimum average number of photons is 
attained for the state 



Vp|ii)-yT^|oo), (14) 

but its entanglement after decoherence is only E « 
0.1622. (One notices the similarity of the state 1)1411 and 
the two-mode squeezed state.) These two facts clearly re- 
fute the intuition that a lower number of photons should 
lead to less decoherence. On the other hand, it is true 



that the state with the largest amount of photons, 4/3, 
decoheres the most. That is, in the state 

yT^|ii) + ^|oo) (15) 

the amount of entanglement surviving the 50% photon 
absorption is £^ ~ 0.0438. The minimum entanglement 
after decoherence for a fixed average number of photons is 
in fact obtained for the family of states for which a — j3. 

In order to explain the preceding results, let us give an 
intuitive idea for why the state 

v/T72|ll) + V172IOO) (16) 

is less robust (when using the entanglement of formation 
as entanglement measure) than the state 

v/T72|10) + v^|01), (17) 

although both states contain 1 photon on average. For 
concreteness, again assume we lose 50% of the photons. 
Suppose we would be able to measure the number of pho- 
tons in the two environments into which the two modes 
decohere (with a perfect photon detector). In that case, 
if we don't find any photons we know we have an en- 
tangled state of the two modes, whereas if we do find a 
photon in the environment, we have no entanglement left 
between the two modes. In the case of the state we 
have a probability of (l-|-77^)/2 = 5/8 to find no photons 
in the environment, but the state of the two modes is 
collapsed onto the (mmormalized) state 

^/i72|ll) + v/T72?y|oo) (18) 

with an entanglement oi E 0.72 for 77 = 1/2. On the 
other hand, for the state (|17|l we have a slightly smaller 
chance of = 1/2 to find no photons in the environment, 
but the state of the two modes then collapses back onto 
()17)l . with its full entanglement of 1 ebit. Thus on aver- 
age we indeed retain more entanglement from state (|17f) 
(namely 0.5 ebits) than from H16() (namely, 0.45 ebits) af- 
ter decoherence, making (I17II the more robust state. This 
then also indicates why the family of states 1)13)1 are more 
robust than either of the families of states 1)14)1 and 1)15)1 . 

A second type of states for which the entanglement 
of formation after decoherence can be calculated is the 
set of entangled coherent states. The reduced density 
matrix of those states has rank two, the support given 
by the two-dimensional Hilbert space spanned by two 
different coherent states, which we choose here as \a) 
and I — a). After decoherence the relevant Hilbert space 
remains two-dimensional, spanned by coherent states | ± 
^/rja). We can distinguish three types of states. The first 
two are both symmetric under the interchange of the two 
modes A and B, 

Vp|+>A|+>i3 + v/T^exp(i</.)|-)A|-)B, (19) 
with < p < 1/2, and 

^ri^\+)A\ + )B + VP^M^4>)\-)A\-)B. (20) 



4 



0.35 




0.5 1 1.5 2 2.5 

number of photons 



FIG. 2: Entanglement of formation as a function of the aver- 
age number of photons for three types of entangled coherent 
states lO, and 

The third type is of the general form 

V^|+>A|->B + yT^exp(z0)|-)A|+>s. (21) 
Here the |±) states are defined by 

|±) = {\a)±\-a))/^ 

N± ^ 2±2exp(-2|an, (22) 

which form an orthonormal basis for the relevant Hilbert 
space. Just as before we choose p = 1/3 and i] = 1/2. 
The entanglement left after decoherence is plotted in 
Fig. 2 for various values of the phase 0. One clearly 
identifies three groups of curves. They correspond to 
the states (PUll . and (|^ . In particular, the states 

with the largest amount of entanglement correspond to 
(|21|l . and the states with the smallest and largest num- 
bers of photons, respectively, to H19() and H20(l . Within 
each group the states have the property that fewer pho- 
tons lead to less decoherence. On the other hand, the 
state with the largest amount of entanglement is not the 
state with the smallest number of photons. The three op- 
timal states within the three groups are the same three 
states singled out above, ((T^ . and ifT^ . and corre- 
pond to the limit a — s- 0. We can infer two conclusions 
from the results plotted in Figs. 1 and 2 that make it seem 
unlikely the two-mode squeezed state would be the most 
robust. First, states with fewer photons are not more ro- 
bust, and second, the symmetric states of the two modes 
are not the most robust either. 

Nevertheless, if we consider the two-mode squeezed 
state with the same amount of entanglement as the other 
states discussed in this subsection, then we find, using 
the formula for the entanglement of symmetric Gaussian 
states from 0, the average number of photons in the 
state and the entanglement left after decoherence are, 
respectively, n « 0.5138 and E 0.3979, thus clearly 
improving on the most robust states considered above. 
Unfortunately, without being able to calculate the en- 
tanglement of formation for more complicated states or 



states with fewer symmetry properties or living in larger 
Hilbert spaces, it is not possible to conclude anything 
yet about the two-mode squeezed state being the most 
robust entangled state of light in this context. That is 
the main reason to consider a measure of entanglement 
that can be calculated easily for any type of entangled 
state, the negativity 

B. Negativity 

Suppose we have an arbitrary pure state of two modes 
A and B. We can always write this state in its Schmidt 
decomposition, 

|^'>AB (23) 

k 

with the Cfc real positive coefficients with ~ 1, and 

\(f>k)A and \ipk) B orthonormal states on systems A and B. 
The negativity of a pure state can be defined in terms of 
the Schmidt coefficients as 

N=\{^c,Y^l). (24) 

k 

For a mixed state the negativity is determined by the 
sum of the absolute values of the negative eigenvalues of 
the partial transpose of the density matrix j^] ■ Numer- 
ically, we expand each density matrix in number states 
and then the partial transpose (with respect to system 
A) is defined as 

{jiA,nB\p^^\mA,mB) = (m^, rislplnA, ms), (25) 

and the negativity is then 

with 1 1 . 1 1 denoting the trace norm. We are interested in 
fixing the entanglement of a pure state and then calcu- 
lating the negativity of the decohered state. Since we 
have now considered two different measures of entangle- 
ment, we will in fact always fix both the entanglement 
of formation and the negativity of the pure states. This 
is done as follows: suppose we are looking for a state 
with M Schmidt coefficients ci . . . cm with M > 4 and 
a fixed value E for the entanglement of formation and a 
fixed value N for the negativity. If we treat C3 . . . Cj\/ as 
fixed, then ci and C2 are determined by the normalization 
and N . Since both the negativity and the norm are sim- 
ple functions of the Schmidt coefficients we can in fact 
determine ci and C2 analytically. The entanglement of 
formation is not a simple function, and so we use a nu- 
merical method to determine coefficients ci . . . C3 given 
coefficients C4 . . . cm- Namely, we use Newton's method 
to obtain better and better estimates of C3, such that the 
entanglement of formation approaches where for each 
C3 the coefficients ci and C2 are fixed analytically. 



0.25 



0.226 




0.224 - 



0.515 0.52 0.525 0.53 

number of photons 

FIG. 3: Negativity as a function of average number of pho- 
tons for 1000 randomly generated states whose Schmidt coef- 
ficients are equal to that of a two-mode squeezed state, with 
a restriction on the distance from a two-mode squeezed state, 
as explained in the text. The circle indicates the two-mode 
squeezed state. 



We first consider the most straightforward case, where 
we in fact keep all Schmidt coefhcients the same as for 
a two-mode squeezed state, but vary the basis vectors 
\(t>k)A and |<^fc)s in the Schmidt decomposition H24|) . An 
arbitrary M-dimensional basis can be fully be specified 
by M — 1 angles, in analogy to the 2 Euler angles in 3-D 
space. We choose the definitions of the angles such that 
the number state basis |0), . . . \ AI — 1) corresponds 
to setting all angles equal to zero. Instead of varying 
those basis vectors over all possible choices by varying 
all angles over 2tt independently, we instead choose, for 
computational efficiency reasons, a large number of ran- 
domly chosen basis vectors by choosing randomly M — I 
angles. In order to choose random states that are within 
a certain "distance" from the standard Fock basis we 
simply restrict the maximum values of the M — 1 angles. 
By increasing the distance one can get further and fur- 
ther away from the two-mode squeezed state. In Figs 3 
and 4 we plot the results. We plot the negativity of the 
decohered state as a function of the average photon num- 
ber in the pure state for two different distances. The first 
restricts the angles to be less than 1/10, in Fig. 4 the dis- 
tance is unrestricted. The two-mode squeezed state we 
take has the same entanglement of formation as before, 
and 1] — 0.5, also as before. The Schmidt coefficients 
of the two-mode squeezed state become exponentially 
smaller with the number of photons n (see Hence, 
if we truncate the Hilbert space at a certain maximum 
photon number iV„i, the state thus obtained is a good ap- 
proximation to a two-mode squeezed state provided Nm 
is sufficiently large. Here we chose N„i = 6. This means 
that by varying the basis states, the maximum possible 
number of photons is, of course, at most 6, as indeed is 
clear in Fig. 4. What we see from the two figures is that 
indeed the two-mode squeezed state itself is the most ro- 
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FIG. 4: Same as previous Figure, except there is no restriction 
on the distance to the two-mode squeezed state. 



bust among the set of states considered here, which can 
be confirmed numerically by decreasing the distance even 
more than in Fig. 3. 

It is probably worthwhile noting that all entanglement 
properties of a given state are determined by the Schmidt 
coefFcients. And so the conclusion is that, with all entan- 
glement properties being equal, the state with the small- 
est number of photons is the most robust, provided we 
use the negativity of the decohered state to measure ro- 
bustness. 

As an aside we note it might be slightly confusing to 
read in Ref. that the negativity of any state, mixed 
or pure, is proportional to the "robustness" of entangle- 
ment. However, the two definitions of the word "robust- 
ness" are quite different. In the robustness of the 
entanglement of a state refers to how much of unentan- 
gled states on the M-dimensional Hlbert space has to 
be mixed with the original state in order to remove all 
entanglement from the state. 

We can try to confirm the behavior of the negativity 
for the same set of states that featured in Fig. 1, i.e., the 
states with at most one photon in each mode A and B. 
One sees that the negativity, in contrast to the entangle- 
ment of formation (See Fig. 1), does seem to favor states 
with fewer photons as far as robustness is concerned. 

On the other hand, there are states that do have the 
same amount of entanglement of formation and the same 
negativity as the two-mode squeeze state, but do not have 
the same Schmidt coefficients. For those states there nec- 
essarily exist at least one measure of entanglement that is 
different than that for the two-mode sqeezed state. Let us 
investigate this a little further. If we restrict the Hilbert 
space of both modes A and B to contain at most three 
photons, then there are at most four Schmidt coefhcients. 
The coefficients are then fixed by the normalization, by 
the negativity, by the entanglement of formation, and 
by some third measure of entanglement, i.e. some other 
independent function of the Schmidt coefhcients Ck- For 
instance, let us, quite arbitrarily, choose the purity of the 
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FIG. 5: Negativity as a function of the average number of 
photons for states of the form @. 
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FIG. 6: Negativity as a function of the average number of 
photons for states with fixed pure-state negativity and entan- 
glement of formation and at most containing three photons 
per mode. 

reduced density matrix as the third entanglement mea- 
sure, 

P = l-Y.i^t). (27) 

We then take a two-mode squeezed state but truncate its 
Fock-state expansion after four terms. In order for it to 
be a good approximation to an actual two-mode squeezed 
state we choose the average number of photons and the 
entanglement smaller than we did in previous examples. 
In Fig. 6 we plot the negativity of the decohered states 
(again using 77 = 0.5) as a function of the average num- 
ber of photons for 1000 randomly generated states that 
all have the same pure-state entanglement E — 0.2 and 
negativity N — 0.2079. The plot then shows that there 
is indeed a line of points, indicating there is exactly one 
more function characterizing the initial pure state. In 
Fig. 7, then, we plot the same data but as a function 
of the purity P, which then indeed shows a monotonic 
behavior, rather than the bistable behavior of Fig. 7. In 
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FIG. 7: Negativity as a function of the purity for states with 
fixed pure-state negativity and entanglement of formation and 
at most containing three photons per mode. 
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FIG. 8; Same as Fig. 6 except for states with at most four 
photons. 

fact, it shows, perhaps surprisingly, that the smaller the 
purity, the more robust the state. From Figs 6 and 7 it 
follows already that the two-mode squeezed state is not 
the most robust entangled state of light, although it is 
close in the sense that the entanglement of the most ro- 
bust state is only a trifle larger than that of the two-mode 
squeezed state. 

We extend this discussion to states confined to rank 
five, containing at most four photons. The correspond- 
ing plots are in Figs. 8 and 9. Those plots show that, 
as expected, there is now more than one extra degree 
of freedom (in this case exactly two degrees of freedom, 
of course) that determines how robust a given entangled 
state is. But the extra degrees of freedom, in this exam- 
ple, hardly increase (less than 0.1%) the negativity of the 
decohered state relative to that of the two- mode squeezed 
state. This is because the states we considered have only 
a small projection onto the extra dimensions (i.e., the 
Fock states with 3 and 4 photons). If we increase the 
amount of entanglement and, concomittantly, the allowed 
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FIG. 9: Same as Fig. 7 except for states with at most four 
photons. 

average number of photons in the initial state, the extra 
degrees of freedom should become more important. For 
example, taking states with an entanglement of E Ri 0.82 
leads to an average number of photons of n « 0.5137 for 
the two-mode squeezed state, and the most robust state 
then turns out to be 2% more robust than the two-mode 
squeezed state. The optimum state has in that case an 
average photon number of n « 0.5797, as was found nu- 
merically. Thus, we expect the two-mode squeezed state 
© to be close to the most robust entangled state of light 
if its entanglement and average number of photons are 
small. For reasonable (i.e., experimentally achievable) 
values the robustness is always clearly within 1% of the 
optimum. 

IV. DISCUSSION 

In this paper we considered the degradation of entan- 
glement of entangled states of light suffering from photon 



absorption losses. The aim was to find the most robust 
entangled state, i.e., the state that preserves its entangle- 
ment best under a given amount of noise. We found that 
the answer depends on how one formulates the problem. 
In particular, it depends on what measure(s) of entan- 
glement one uses. That is, we may fix any number n 
of entanglement measures Mi . . . M„ of the initial pure 
state and use one particular entanglement measure Mq of 
the decohered state (after photon absorption) to define 
the robustness of the state. The answer which state is 
the most robust then depends on the measures Mq and 
Ml . . . Mn one uses. 



If one fixes two measures of entanglement and takes Mi 
to be the entanglement of formation and M2 the nega- 
tivity, and one uses the negativity for Mq, we found that 
states with smaller number of photons tend to be more 
robust. Nevertheless, the most robust state is not the 
state with the smallest possible number of photons, the 
two-mode squeezed state. The reason is that even when 
one fixes both the entanglement of formation and the 
negativity there are still other degrees of freedom left 
that determine the entanglement properties of the pure 
state, including its robustness against photon absorption 
losses. The two-mode squeezed state is close to the most 
robust state, though, and is very close for experimentally 
achievable parameters. That is, the value of Mq for the 
most robust state differs less than 1% from that of the 
two-mode squeezed state for realistic parameters. 



On the other hand, if one fixes all entanglement degrees 
of freedom of a pure state (i.e. all Schmidt coefficients) 
and uses the negativity of the decohered state to quantify 
robustness, then the most robust state is the one with the 
smallest (given all constraints) number of photons. 
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